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Brownian motion has played important roles in many different fields of science since its origin 
was first explained by Albert Einstein in 1905. Einstein's theory of Brownian motion, however, is 
only applicable at long time scales. At short time scales, Brownian motion of a suspended particle 
is not completely random, due to the inertia of the particle and the surrounding fluid. Moreover, 
the thermal force exerted on a particle suspended in a liquid is not a white noise, but is colored. 
Recent experimental developments in optical trapping and detection have made this new regime 
of Brownian motion accessible. This review summarizes related theories and recent experiments 
on Brownian motion at short time scales, with a focus on the measurement of the instantaneous 
velocity of a Brownian particle in a gas and the observation of the transition from ballistic to diffusive 
Brownian motion in a liquid. 



I. INTRODUCTION 

Brownian motion is the apparently perpetual and ran- 
dom movement of particles suspended in a fluid (liquid or 
gas), which was first observed systematically by Robert 
Brown in 1827 When Brown used a simple micro- 
scope to study the action of particles from pollen im- 
mersed in water [1], he "observed many of them very 
evidently in motion". The size of those particles was 
about 5 /im. He also observed the same kind of motion 
with powders of many other materials, such as wood and 
nickel, suspended in water. 

As first explained by Einstein in 1905 [2], the Brow- 
nian motion of a suspended particle is a consequence of 
the thermal motion of surrounding fluid molecules. Ein- 
stein's theory of Brownian motion predicts that 

([Ax(t)} 2 ) = ((x(t) - (x(0)) 2 ) = 2Dt, (1) 

where ([Ax(t)] 2 } is the mean-square displacement (MSD) 
of a free Brownian particle in one dimension during time 
t , and D is the diffusion constant. The diffusion constant 
can be calculated by D = ksT/j, where ks is the Boltz- 
mann constant, T is the temperature, and 7 = 6ttt]R is 
the Stokes friction coefficient for a sphere with radius R. 
Here 77 is the viscosity of the fluid. 

M. von Smoluchowski also derived the expression of 
MSD independently in 1906 [3[, with a result that dif- 
fered from Eq. (pQ) by a factor of about 2. In 1908, Paul 
Langevin introduced a stochastic force and derived Eq. 
(PQ) from Newton's second law Langevin's approach 

is much more intuitive than Einstein's approach, and the 
resulting "Langevin equation" has found broad applica- 
tions in stochastic physics [6[ . Experimental confirmation 
of Eq. (pQ) was provided by the brilliant experiments of 
Jean Perrin [7[ , recognized by the Nobel Prize in Physics 
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in 1926. Theodor Svedberg also verified the Einstein- 
Smoluchowski theory of Brownian motion and won the 
Nobel Prize in Chemistry in 1926 for related work on 
colloidal systems [8[. 

Persistence and randomness are generally accepted as 
two key characteristics of Brownian motion. The tra- 
jectories of Brownian particles are classic examples of 
fractals [9[. They are commonly assumed to be con- 
tinuous everywhere but not differentiable anywhere fioj ] . 
Since its trajectory is not differentiable, the velocity of 
a Brownian particle is undefined. According to Eq. (pp), 
the mean velocity m easured over an interval of time t 
is v = yJ([Ax(t)] 2 )/t = V2D/VI. This diver ges as t 
approaches 0, and therefore does not represent the real 
velocity of the particle [ll|, [l2| . 

At short time scales (t <C r p , where r p = M/7 is the 
momentum relaxation time of a particle with mass M), 
the dynamics of a Brownian particle is expected to be 
dominated by its inertia and its trajectory cannot be self- 
similar. This is termed "ballistic Brownian motion" to 
be distinguished from the common "diffusive Brownian 
motion" . Fig. [1] shows a 2D trajectory of a Brownian 
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FIG. 1. A 2D trajectory of a Brownian particle. The black 
curve is assumed to be a true trajectory of the particle. Red 
dots are measured positions, and red curves are measured 
trajectories. The sampling rate of (B) is 10 times of that of 
(A). 
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particle. The black curve is assumed to be a true trajec- 
tory of the particle. Red dots are measured positions. In 
Fig. [T]A., the sampling rate is too small to measure the ve- 
locity of the Brownian particle. The measured trajectory 
(red curve) is completely different from the true trajec- 
tory, and appears chaotic. It is impossible to obtain the 
velocity of the particle from the measured trajectory in 
Fig. [TJA.. In Fig. [TJ3, the sampling rate is much larger. 
Now the measured trajectory is very close to the true 
trajectory of the particle. If the measured displacement 
of the particle is Ax(t) during time At, then the velocity 
of the particle is approximately v = Ax(t)/At. 

In 1900, F. M. Exner made the first quantitative study 
of Brownian motion by measurin g th e velocity of Brown- 
ian particles suspended in water [13l [l4j . He found that 
the measured velocity decreased with increasing parti- 
cle size and increased with increasing water temperature. 
However, his measured velocities were almost 1000-fold 
smaller than those predicted by the energy equipartition 
theorem [13j . The reason of this discrepancy was not un- 
derstood until A. Einstein developed his kinetic theory 
about Brownian motion |2j. 

In 1907, Einstein published a paper entitled "Theo- 
retical observations on the Brownian motion" in which 
he considered the instantaneous velocity of a Brownian 
particle [ll|, [l2| . Einstein showed that by measuring this 
quantity, one could prove that "the kinetic energy of the 
motion of the centre of gravity of a particle is independent 
of the size and nature of the particle and independent of 
the nature of its environment" . This is one of the basic 
tenets of statistical mechanics, known as the equiparti- 
tion theorem. However, Einstein concluded that due to 
the very rapid randomization of the motion, the instanta- 
neous velocity of a Brownian particle would be impossible 
to measure in practice [ll|, [l2| : 

"We must conclude that the velocity and direction of 
motion of the particle will be already very greatly altered 
in the extraordinary short time and, indeed, in a to- 
tally irregular manner. It is therefore impossible - at 

least for ultramicroscopic particles - to ascertain \frP by 
observation." 

Einstein's conclusion was unchallenged for more than 
100 years because of the exclusive difficulty of such a 
measurement. For a 1 jam diameter silica (SiCte) sphere 
in water at room temperature, the momentum relaxation 
time t p is about 0.1 /is and the root mean square (rms) 
velocity v rms = y^ksT/M is about 2 mm/s in one dimen- 
sion. To measure the instantaneous velocity with 10% 
uncertainty, one would require 2 pm spatial resolution 
and 10 ns temporal resolution, which is a very difficult 
task. Due to the lower viscosity of gas as compared to 
liquid, the momentum relaxation time r p of a particle in 
air is much larger. This lowers the requirements of both 
temporal and spatial resolution. 

Nondiffusive Brownian motion of colloidal suspensions 
with high concentrations at short time scales have been 
studied by measuring the autocorrelation functions of 
multiply scattered, transmitted light [l5l-[l7|. Recent 



developments in optical tweezers and detection systems 
with unprecedented resolution now prove to be an in- 
dispensable tool for studying the Brownian motion of a 
single particle at short time scales [l8l-[26j. For exam- 
ple, we have observed the ballistic Brownian motion and 
measured the instantaneous velocity of a Brownian par- 
ticle for the first time with an optically trapped bead in 
air [2l|. Huang et al have observed the transition from 
ballistic Brownian motion to diffusive Brownian motion 
in a liquid [24| . Franosch et al have observed resonances 
arising from hydrodynamic memory in Brownian motion 
in a liquid and the long-sought colored spectrum of the 
thermal force [H,^. 

Besides Brownian motion at short time scales, theo- 
retical and experimental studies of anisotropic Brownian 
motion and Brownian motion in nonequilibrium systems 
are currently pursued by many groups. For example, 
several groups reported anisotropic Brownian motion of 
particles near interfaces [28-32], and Brownian motion 
of anisotropic particles such as ellipsoids [33|, nan- 
otubes [35|, [36| and helical bacteria [37j • Brownian mo- 
tion in nonequilibrium systems is of particular interest 
because it is directly related to the transport of molecules 
and cells in biological systems. Important examples in- 
clude Brownian motors [38, 391, active Brownian motion 
of self-propelled particles [40-46], hot Brownian motion 
(47(1 , and Brownian motion in shear flows [48]. Recent 
theoretical studies also found that the inertias of par- 
ticles and surrounding fluids can significantly affect the 
Brownian motion in nonequilibrium systems [49l-l54j . 

In this review, section 2 introduces the theories of 
Brownian motion of particles in a gas, and the recent 
measurement of the instantaneous velocity of a Brown- 
ian particle in air. Section 3 introduces the theories of 
Brownian motion of particles in a liquid at short time 
scales, the experimental observation of the colored ther- 
mal force, and the transition from ballistic to diffusive 
Brownian motion in a liquid. Section 4 discusses the ef- 
fects of detection noise on the measurement of different 
quantities of Brownian motion. Finally, in section 5, we 
discuss future experiments on Brownian motion at short 
time scales. 



II. BROWNIAN MOTION IN A GAS 

In this section, we assume that the density of a gas is 
much smaller than the density of the suspended Brow- 
nian particles. So the inertia effects of the gas can be 
neglected. 



A. Theory 

The mean free path of molecules in air at 1 atmosphere 
at room temperature is about 68 nm [55|. The collision 
rate between a l-/im- diameter microsphere suspended in 
air and surrounding air molecules is about 10 16 Hz at 
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ambient conditions. The observed Brownian motion is 
an averaged effect of these ultrafast collisions. Because 
of the huge difference between the mass of a microsphere 
and that of an air molecule, the motion of a microsphere 
can only be changed significantly by a large number of 
collisions. This is reflected in the fact that r p = 6 /is 
for a l-/im-diameter silica microsphere in air at ambient 
conditions. The dynamics of a Brownian particle at time 
scales much longer than that of individual collisions can 
be described by a Langevin equation [4, 5]. Here we in- 
troduce the theory of Brownian motion in a gas following 
the classic work of Uhlenbeck and Ornstein [56j . 



1. A free particle in a gas 

The dynamics of a Brownian particle with mass M in 
a gas can be described by a Langevin equation 0, E| EH • 



where rj is the viscosity coefficient of the gas, R is the 
radius of the microsphere, and Kn = s/R is the Knud- 
sen number. Here s is the mean free path of the gas 
molecules. c K = (0.31Kn)/(0.785 + 1.152Kn + Kn 2 ) is 
a small positive function of Kn [58[. At high pressures 
where KnC 1, the damping coefficient is To = 67rr]R/M, 
which is the same as the prediction of Stokes' law. 

At very short time scales, the motion is ballistic and its 
instantaneous velocity can be measured as v = Ax(i)/t, 
when t < r p [56]. The ballistic Brownian motion is dif- 
ferent from a simple ballistic motion. For a simple bal- 
listic motion with velocity u, we have Ax(t) = ut and 
[Ax(t)] 2 = u 2 t 2 . The velocity u can be any value and 
usually has no relation with the temperature of the envi- 
ronment. For the ballistic Brownian motion, the ampli- 
tude of the velocity is determined by the temperature of 
the environment. The ID Maxwell-Boltzmann distribu- 
tion of the velocity of a particle in thermal equilibrium 



where 



F therm (t) = (2fc B T 7 ) 1 / 2 C(t) 



(2) 



(3) 



is the Brownian stochastic force. is a normalized 

white- noise process. Hence for all t and t' ', 



(C(t)> = , and (at)C(t')) = 5(t-t'). 



(4) 



The MSD for a Brownian particle at thermal equilib- 
rium with the air is [56]: 



2k B T 
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where Tq = j/M is the damping coefficient. We have 
t p = 1/Tq. At long time scales (t ^> r p ), the MSD is the 
same as that predicted by Einstein's theory (Eq. (pQ)). 
At very short time scales, the MSD is 



M r for t « r p . 



The velocity autocorrelation function is [561 ]: 
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Although these equations are initially derived for an 
ensemble of particles, the ergodic theorem dictates that 
they are also valid for measurements of a single particle 
taken over a long time. 

The damping coefficient To in a gas can be calcu- 
lated by kinetic theory. Assuming the reflection of gas 
molecules from the surface of a microsphere is diffusive, 
and the molecules thermalize with the surface during col- 
lisions, we obtain [53, l5q 

67rr]R 0.619 M 
ro = ^0619T^ (1+C ^ (8) 
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where V{ is the velocity of the particle along direction z, 
which can be any direction. 



2. An optically trapped microsphere in a gas 

For small displacements, the effect of optical tweez- 
ers on the microsphere's motion can be approximated by 
that of a harmonic potential. The Brownian motion of a 
particle in a harmonic trap has been studied by Uhlen- 
beck and Ornstein [56|, and Wang and Uhlenbeck (59^ . 
The equation of the Brownian motion of a microsphere 
in a harmonic trap is: 



d x 



dx 
] ~dt 



tt 2 x 



(10) 



where £1 = y / 1 Kjm is the natural angular frequency 
of the trapped microsphere when there is no damping, 
and A = (2k B TT /M) 1 ^ 2 . The cyclic frequency of the 
damped oscillator is uj\ = V^ 2 - r§/4. The system 
is underdamped when uj\ is real (Q > T/2), critically 
damped when uji = 0, and overdamped when oj\ is imag- 
inary (fi < T/2). 

The MSD of a Brownian particle in an underdamped 
harmonic trap in air is |59j : 



([Ax(t)] 2 ) = 



2k B T 
Mft 2 



-t/2r p 



COSUJit 



sinc^it 

2uJ\T p 



(11) 

The position autocorrelation function is related to the 
MSD by: 

([Ax(t)} 2 )=2(x 2 )-2(x(t)x(0)), (12) 

where (x 2 ) = k B T '/ \MVt 2 ). The rms amplitude is x rms = 
^k B T I (MQ 2 ). The normalized position autocorrelation 
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function (PACF) of the particle is [59|: 



1. A fast detection system 



(x(t)x(0)} 



= e-*/ 2T 



COSUJit - 



sinc^it 

2(jJiT p 



(13) 



The normalized velocity autocorrelation function 
(VACF) of the particle is 



(v(t)v(0)) 
(v 2 ) 



-t/2r p 



COSUJit 



sinc^it 



(14) 



Both the position autocorrelation function and the veloc- 
ity autocorrelation function oscillate for an underdamped 
system. 

Similar to the optical spectrum of an atom, the power 
spectrum of the Brownian motion of a trapped micro- 
sphere contains a lot of information about the system. 
The power spectral density (PSD) of a variable is the 
squared modulus of its Fourier transform (59l46lj. The 
expected values of the PSD of an optically trapped mi- 
crosphere in air is 



2k B T 



Mn 2 (n 2 - w 2 ) 2 + w 2 r 2 



2 • 



(15) 



Eq. ([T5]) is valid for both underdamped [57j and over- 
damped systems. The measured PSD of a recorded 
x(t) is 5^ ec , which depends on a sample of the white 
noise ((t). Thus an experimental PSD will appear noisy. 
Averaging many measured 5^ ec will result in a spec- 
trum close to the expected spectrum S(uj). Another 
way to reduce the noise in a measured spectrum is 
"blocking" [60] . A "block" of consecutive data points 
(cj^ ,5^ c )... {^k 2 ^kl c ) can be represented by a single 
new "data point" (uJk^S 1 ^ ) which are the block averages. 



In order to measure the instantaneous velocity of the 
Brownian motion of a trapped microsphere in air, we 
implemented an ultrahigh resolution detection system. 

We used an ultra-stable NPRO laser (Model: 126- 
1063-700, Lightwave Electronics (now JDSU)) to trap 
and monitor a silica bead in vacuum. Its rms inten- 
sity noise is < 0.05 % over the range from 10 Hz to 2 
MHz, and is shot noise limited above 10 MHz. It is a 
single frequency laser with a wavelength of 1064 nm and 
a coherence length longer than 1000 m. A detailed char- 
acterization of this type of laser can be found in Ref. [64j . 
We used this laser for both trapping and detection. This 
is achieved by using a polarizing beam splitter cube to 
reflect one of the trapping beams for detection (Fig. [2j). 

Our lab has previously developed a fast position- 
sensitive laser beam detector [20]. The previous detec- 
tor used a bundle of optical fibers that spatially splits 
the incident beam, and a fast balanced photodetector 
to measure the difference between the two halves of the 
beam. We simplified the detection system by using a mir- 
ror with a sharp edge (BBD05-E03, Thorlabs) to replace 
the fiber-optic bundle for splitting the beam (Fig. [2]). 
The sharp edge of a mirror is much smoother than the 
boundary between the two halves of a fiber bundle. So 
it is much simpler and has less noise than a fiber bundle 
for splitting the laser beam. 

We used a balanced detector (PDB120C, photodiode 
diameter: 0.3 mm, Thorlabs) with a bandwidth of 75 
MHz for detection. The detector is sensitive to light 
with wavelengths in the range of 800-1700 nm. It has 
a high transimpedance gain of 1.8 x 10 5 V/A. The de- 
tector measures the difference between the two halves of 



B. Experimental observation of the instantaneous 
velocity of a Brownian particle in air 

Because of the lower viscosity of air as compared to 
that of liquid, a particle suspended in air is an ideal 
system for studying the ballistic Brownian motion. The 
main difficulty of studying Brownian motion in air, how- 
ever, is that the particle will fall under the influence of 
gravity. To overcome this problem, Fedele et al used 
aerosol suspensions with small particles (~ 0.2 fim diam- 
eter) to achieve a long sedimentation time [62|, Blum et 
al performed the experiment under microgravity condi- 
tions in the Bremen drop tower [63[ . Due to the lack of a 
detection system with a sufficient resolution, neither ex- 
periment was able to measure the instantaneous velocity 
of Brownian motion. We overcame this problem by us- 
ing optical tweezers to simultaneously trap and monitor 
a silica bead in air and vacuum, allowing long-duration, 
ultra-high-resolution measurements of its motion. Here 
we review this experiment that was originally reported 
in Ref. [2l| with more experimental details. 



Vacuum 
Chamber 




FIG. 2. This simplified schematic shows our counter- 
propagating dual-beam optical tweezers, and a novel detec- 
tion system (Figure adapted from Ref. [2lj]). The s-polarized 
beam is reflected by a polarizing beam splitter cube after 
it passes through a trapped bead inside a vacuum chamber. 
Then, for detection, it is split by a mirror with a sharp edge. 
The ^-polarized beam passes through the cube. 
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the beam, which is proportional to the particle excursion. 
The intensity noise of the laser is contained in both halves 
and is thus canceled in the measurement. This detection 
system enables us to monitor the real-time position of a 
trapped microsphere in air with Angstrom spatial resolu- 
tion and microsecond temporal resolution [2l| . We have 
also developed a 3D detection system that can monitor 
the 3D motion of a microsphere trapped in vacuum with 
a sensitivity of about 39 fm Hz _1//2 over a wide frequency 
range [57j. 



2. Experimental results 

A simplified scheme of our setup for measuring the 
instantaneous velocity of a Brownian particle in air is 
shown in Fig. [2] The trap is formed inside a vacuum 
chamber by two counter-propagating laser beams focused 
to the same point by two identical aspheric lenses with 
focal length of 3.1 mm. The two 1064 nm laser beams 
are orthogonally polarized, and their frequencies differ by 
160 MHz to avoid interference. 

The two laser beams are aligned with the help of a 
pinhole aperture whose diameter is 1.0 ± 0.5 jam. We 
intentionally make the waist of one beam larger than the 
other to make this alignment less critical. The measured 
waists of the two beams are about 2.2 jam. and 3.0 /im, 
respectively. Once a bead is trapped, we keep the power 
of one beam constant, and tune the power of the other 
beam to maximize the trapping frequency. 

When the bead deviates from the center of the trap, 
it deflects both trapping beams. We monitor the posi- 
tion of the bead by measuring the deflection of one of the 



beams, which is split by a mirror with a sharp edge. This 
simple, yet novel, detection scheme has a bandwidth of 
75 MHz and ultra- low noise [20|, [65|. The position sig- 
nal of a trapped bead is recorded at a sampling rate of 2 
MHz. Because of the detection noise, we are not able to 
obtain accurate instantaneous velocities of a bead at this 
rate. To reduce the noise, we average every 10 successive 
position measurements, and use these averages to calcu- 
late instantaneous velocities with time resolution of 5 /is. 
Although this method reduces the temporal resolution 
by a factor of 10, it greatly increases the signal-to- noise 
ratio if both the trapping period (2tt/ujo) and momen- 
tum relaxation time are much larger than 5 /is. These 
conditions are satisfied here since the trapping period is 
about 320 /is, r p = 48 /is at 99.8 kPa (749 torr), and 
r p = 147 fjs at 2.75 kPa (20.6 torr). 

Figure [3] shows typical samples of position and velocity 
traces of a trapped bead. The position traces of the bead 
at these two pressures appear very similar. The instan- 
taneous velocity of the bead at 99.8 kPa changes more 
frequently than that at 2.75 kPa, because the momentum 
relaxation time is shorter at higher pressure. 

Figure [4] shows the mean square displacements of a 
3 fim silica bead as a function of time. The measured 
MSD's fit excellently with Eq. [TTIover three decades of 
time for both pressures. The measured MSD's are com- 
pletely different from those predicted by Einstein's theory 
of Brownian motion in a diffusive regime. The slopes of 
measured MSD curves at short time scales are double of 
those of the MSD curves of diffusive Brownian motion 
in the log- log plot (Fig. [4]). This is because the MSD is 
proportional to t 2 for ballistic Brownian motion, and it is 
proportional to t for diffusive Brownian motion. Another 
important feature is that the MSD curves are indepen- 
dent of air pressure at short time scales, as is predicted 
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FIG. 3. One-dimensional trajectories of a 3 fim diameter silica 
bead trapped in air at 99.8 kPa (A) and at 2.75 kPa (B). The 
instantaneous velocities of the bead corresponding to these 
trajectories are shown in (C) and (D). Figure adapted from 
Ref. El. 
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FIG. 4. The mean square displacements of a 3 fim silica 
bead trapped in air at 99.8 kPa (red square) and 2.75 kPa 
(black circle) . The solid lines are the theoretical predictions of 
Eq. Illl The prediction of Einstein's theory of free Brownian 
motion in the diffusive regime is shown in dashed lines for 
comparison. Figure adapted from Ref. [2 II ]. 
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by Eq. [6]for ballistic Brownian motion, whereas the MSD 
in the diffusive regime does depend on the air pressure. 
At long time scales, the MSD oscillates and saturates at 
a constant value because of the optical trap. 

The distributions of the measured instantaneous ve- 
locities are displayed in Fig. El They agree very well 
with the Maxwell-Boltzmann distribution. The mea- 
sured rms velocities are v rms = 0.422 mm/s at 99.8 kPa 
and v rms = 0.425 mm/s at 2.75 kPa. These are very 
close to the pre diction o f the energy equipartition the- 
orem, v rms = s^ksT/M , which is 0.429 mm/s. As ex- 
pected, the velocity distribution is independent of pres- 
sure. Thus the Maxwell-Boltzmann distribution of veloc- 
ities and the equipartition theorem of energy for Brown- 
ian motion were verified directly. 

Figure [6] shows the normalized VACF of the bead at 
two different pressures. They fit with Eq. [14] nicely. At 
2.75 kPa, one can clearly see the oscillations due to the 
optical trap. Eq. [14]is independent of the calibration fac- 
tor of the detection system. The only independent vari- 
able is time t, which we can measure with high precision. 
Thus the normalized VACF provides an accurate method 
to measure r p and c^o • We can also calculate the diameter 
of the silica bead from the r p value at 99.8 kPa [66]. The 
obtained diameter for this microsphere is 2.79 jam. This 
is within the uncertainty range given by the supplier of 
3.0 /am silica beads. We use this value in the calculation 
of MSD and normalized VACF. 

For a particle at a certain pressure and temperature, 
t p should be independent of the trapping frequency. We 
verified this by changing the total power of the two laser 
beams from 25 mW to 220mW. The measured r p of a 
microsphere trapped at 19.6 torr and 749 torr as a func- 
tion of the total laser power is shown in Fig. [71 Here we 
use a new microsphere and a smaller data set to calculate 
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FIG. 6. The normalized velocity autocorrelation functions of 
the 3 fim bead at 99.8 kPa (red square) and at 2.75 kPa (black 
circle) from the measurements. The solid lines are fittings 
with Eq. [TU Figure adapted from Ref. [2lj ] . 
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FIG. 5. The distribution of the measured instantaneous veloc- 
ities of a 3 /im silica bead. The statistics at each pressure are 
calculated from 4 million instantaneous velocities. The solid 
lines are Maxwell-Boltzmann distributions. Figure adapted 
from Ref. 



FIG. 7. Measured momentum relaxation times (r p ) of a mi- 
crosphere trapped at 2.61 kPa and 99.8 kPa as a function of 
the total power of the two trapping beams. 



the r p than those used in the previous figures. Although 
the data points for each pressure are not perfectly on 
a line, it is clear that the r p 's are independent of the 
laser power within the experimental uncertainty. Fitting 
the data for each pressure with a straight line, we ob- 
tain r p = [151.3 + 0.00168 P/(lmW)] /is at 2.61 kPa, and 
r p = [53.74 - 0.00275 P/(lmW)] /is at 99.8 kPa for this 
microsphere, where P is the total power of the two trap- 
ping beams. Thus r p changes less than 1.3% for both 
pressures when the total laser power is changed from 
to 230 mW. This proves that the fitting method is very 
accurate, and the heating due to the laser beams (which 
would change the viscosity [67, 68] and affect r p ) is neg- 
ligible at these pressures. 
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III. BROWNIAN MOTION IN A LIQUID 

The main difference between the Brownain motion in 
a liquid and that in a gas is the hydrodynamic effects 
of the liquid (46j . The Brownian motion of colloidal 
particles in a liquid at high concentrations have been 
studied with diffusing wave spectroscopy, which requires 
each photon to be scattered many times before reaching 
the detector pJ5l - [l7j . The recent developments in opti- 
cal tweezers provide a new tool for studying the Brown- 
ian motion of single particles with unprecedented preci- 
sion [HI, [20|, HH HH Ell • Meanwhile, precise calibrations 
(force, position, etc.) of optical tweezers demand better 
understanding of the Brownian motion of trapped parti- 
cles [H Mmk • 



A. Theory 

Besides the inertia of the particle itself, the inertia of 
the surrounding liquid is also important for Brownian 
motion of particles in a liquid. The motion of a parti- 
cle will cause long-lived vortices in the liquid that will 
affect the motion of the particle itself. This is the hy- 
drodynamic memory effect of the liquid, which domi- 
nates the dynamics of the particle at short time scales. 
These hydrodynamic memory effects were first studied by 
Vladimir sky in 1945 [71|. In 1960s, several authors found 
in computer simulations that the velocity autocorrelation 
function (VACF) of fluid molecules had a power-law tail 
in the form of t~ 3 l 2 [72|, [zsj, in contrast to the exponen- 
tial decay in a dilute gas. Hinch obtained an analytical 
solution of the VACF for free particles from the original 
Langevin analysis [74] . Clercx and Schram calculated the 
MSD and VACF of a Brownian particle in a harmonic po- 
tential in an incompressible liquid [75| , which can be used 
to describe the Brownian motion of an optically trapped 
microsphere in a liquid directly 
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1. A free particle in a liquid 

The effective mass of the microsphere in an incom- 
pressible liquid is the sum of the mass of the microsphere 
and half of the mass of the displaced liquid [77], [78| : 



M* = M n + -M 



(16) 



where M p = (4/3)ttR 3 p p is the mass of the microsphere, 
Mf = (4/3)7rR 3 pf is the mass of displaced liquid, p p is 
the density of the microsphere, and pf is the density of 
liquid. The energy equipartition theorem needs to be 
modified to: 



-M*(v 2 ) 



I. 



:k B T 



(17) 



where v is the velocity of the microsphere in one di men- 
sion. Thus the rms velocity is v rrns = \/k B T '/M* . Be- 



Diameter 


T p 


r f 


T c 


k 




(/im) 


(ps) 


(ps) 


(ns) 


(/iN/m) 


(lis) 


1.0 


0.11 


0.25 


0.34 


100 


94 


3.0 


1.0 


2.2 


1.01 


33.3 


851 


4.7 


2.45 


5.51 


1.58 


21.3 


2083 


10 


11.1 


25.0 


3.4 


10 


9443 



TABLE I. Characteristic time scales of an optically trapped 
silica microsphere in water at 20 °C. Some examples of the 
spring constant of the optical trap (k) are shown in the 5th 
column. It is assumed to be inversely proportional to the 
diameter of the microsphere when the laser power is constant. 



cause of the memory effect of liquid, the velocity auto- 
correlation function (VACF) of a free particle in a liq- 
uid will not be (v(t)v(0)) = ^^e~ l l Tv as in air, but 



M*MQ)) 

k B T/M* 



a + e Q! + t erfc(a + V^) 



a_e a -^erfc(a_ y/i) 



where 



a± = — 



3±(5-36r p /r / ) 1 / 2 

(1 + 97-p/T/) 



_V2 



(18) 



(19) 



r p = M p /(67rr]R) = ^R 2 p p /r] is the momentum re- 
laxation time of the particle due to its own inertia, 
Tf = R 2 pf/r] characterizes the effect of liquid. Here 77 
is the viscosity of liquid and R is the radius of the micro- 
sphere. 

At long time scales, Eq. fT8l approaches 



(v(t)v(0)) 



1 



for t —> 00. 



k B T/M* t 3 / 2 
At short time scales, Eq. [TBI approaches 
(v(t)v(0)) 



k B T/M* 



exp 



for t -)■ 0, 



(20) 



(21) 



where 
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V5F(1 + 2p p / Pf ) ' 

For a silica microsphere in water, b = 2.03. The normal- 
ized VAC F ap proaches 1 at short time scales as 
ex.~p(—by/t/rf), rather than exp(— t/r p ). Thus the dy- 
namics of the particle is dominated by the hydrodynamic 
effects of the liquid. This is very different from the case 
in air. 



2. An optically 



microsphere in a liquid 



The optical trap provides a harmonic force F t 



trap 



—kx on the microsphere when the displacement of the 
microsphere is small, k = M p Q 2 where Q is the nat- 
ural angular frequency of the trap. Clercx and Schram 



8 



[75j gave analytical solutions for the MSD and VACF of a 
trapped Brownian particle in a liquid, and Berg-S0rensen 
and Flyvbjerg [60| gave a solution for the power spec- 
trum density (PSD) of a trapped Brownian particle in a 
liquid. This section introduces their analytical solutions 
and provides some numerical results to visualize those 
solutions. 

Because the velocity of the Brownian motion of a mi- 
crosphere in liquid is much smaller than the speed of 
sound in the liquid, the fluid motion can be described 
by the linearized incompressible time-dependent Navier- 
Stokes equation. The Langevin equation of the motion 
of a trapped microsphere in an incompressible liquid is 
III: 



M*x(t) = -kx(t) - 6in]Rx(t) 



(22) 



-QR 2 ^Wpfq / (t - t')- yi x{t')dt' + F therm (t). 

J — oo 

The first term after the equal sign of Eq. [22] is the har- 
monic force, the second term is the ordinary Stokes's fric- 
tion, the third term is a memory term associated with the 
hydrodynamic retardation effects of the liquid, and the 
last term is the Brownian stochastic force. 

By the famous fluctuation-dissipation theorem, the 
thermal force is directly related to the frictional force. So 
the hydrodynamic memory of the liquid will affect both 
the thermal force and the frictional force. The thermal 
force is not a white noise, but becomes colored. The 
correlation in the thermal force is HE H3, EE Izli 



(Ftherm(t)Ftherm(O)) 



- 7fcsr vi^" 3/2 ' (23) 



which is very different from a delta function (Eq. (TjTJ)). 
The mean-square displacement of a trapped micro- 



sphere in a liquid is [75|, I8C 
2k B T 



{[Ax(t)Y)trap 



(24) 



k 

2k B T 

M* ziOi - z 2 )(z 1 - z 3 )(z 1 - z 4 ) 
e z ^ erfc(z 2 Vi) 



Z 2 (z 2 


~ Zi)(z 2 ~ Z 3 )(Z 2 


- z 4 ) 




e z\t erfc(23V^) 




Z 3 (z 3 


- zi)(z 3 - z 2 ){z 3 


- z A ) 




e 2; 4 £ erfc(z 4 V^) 





z A (z A - zi)(z A - z 2 )(z A - z 3 ) 

The coefficients Z\, z 2l z 3: and z^ are the four roots of 
the equation [80| 



t p + lr f ) z A - ^TJz 3 + z 2 + - = 0, (25) 



where = 6irr]R/k. For t —> oo, Eq. [Ml approaches 

2k B T 



([AX(OO)] } t rap 



k ' 



The normalized VACF of a trapped microsphere in a 
liquid is |zS E3 



A(t) 



(v(t)v(0)} _ z\ e*?* erfc(^iVt) 
k B T/M* ~ (z! - z 2 )( Zl - Z 3 )(Z! - z 4 ) 

z\ e z ^ erfc(z 2 ^/i) 



(26) 



(Z 4 - Zi)(Z4, ~ Z 2 )(Z4 ~ Z 3 ) ' 



(*2 


- zi)(z 2 - z 3 )(z 2 - 


za) 




z\ e 2; 3 £ erfc(z 3 v / t) 






- z x \z 3 - z 2 )(z 3 - 


za) 




z\ e z ^ erfc(z 4 V^) 
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FIG. 8. Calculated mean square displacement of an optically 
trapped silica microsphere in water at 20 °C. Parameters are 
the same as those in Table [T] 




Frequency (Hz) 

FIG. 9. Calculated power spectra of an optically trapped 
silica microsphere in water at 20 °C. Parameters are the same 
as those in Table [I] 
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S(f) 



The power spectral density is |6(J, I8C 
D 



where 



2tt 2 / 2 



(27) 



\[]JWf 



{4>k/f - ^fJJW^ - f/4> P - f/W f ) 2 + (i + ^7/207 ) 2 

where / is the observation frequency, <pk = l/(2irTk) is 
the corner frequency of the power spectrum due to the 
trap, and (f> p j = l/(27rr Pj /). For / — >> 0, Eq. [23 ap- 
proaches 



5(0) 



2fe B T 7 
& 2 ' 



where 7 = 6ttt]R. 

At t 0, Eq. [26] predicts (v(0)i;(0)> = feT/M*, 
which is different from the energy equipartition theorem 
(v(O)v(O)) = ksT/Mp. This conflict is caused by the as- 
sumption in Eq. [22] that the liquid is incompressible. For 
t < t c , we need to consider the liquid to be compressible. 
Here t c = R/c is the time required for a sound wave to 
travel a sphere radius, where c is the speed of sound in 
the liquid. The effects of compressibility have been stud- 
ied by Zwanzig and Bixon [77]. The normalized velocity 



autocorrelation function at t > 

M*M0)> 



t c is 



A(t) 



1 



k B T/M* 
M f 



(28) 
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FIG. 10. Calculated normalized velocity autocorrelation func- 
tion of an optically trapped silica microsphere in water at 20 
°C. Parameters are the same as those in Table [I] The thin 
solid lines (t < 10 ~ 8 s) are calculated from Eq. 1281 and the 
thick solid lines are calculated from Eq. 1261 The dashed lines 
are exponential decays with r p = 1.0/is, corresponding to a 
3.0 jim microsphere. 
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(29) 
(30) 



At very short time scales t < t c , Eq. ([28]) approaches 
A(0) = 1 + jm~' The short time limit A(0) ^ 1 because 
the normalization factor is hsT/M* in Eq. ([28]) . rather 
than k B T/M p . 

The calculated MSD's of microspheres with different 
diameters in water are shown in Fig. [8] and the cor- 
responding power spectra are shown in Fig. [9] Fig. [10] 
displays the normalized velocity autocorrelation function 
(A(t)) of an optically trapped silica microsphere in water 
at 20 °C. The thick solid lines are calculated from Eq. [26] 
which treats the water as an incompressible fluid. The 
thin solid lines at short time scales (t < 10 -8 s) are cal- 
culated from Eq. [28j which includes the compressibility 
effects of water. The dashed lines are exponential decays 
with r p = 1.0 /is, corresponding to a 3.0 fim microsphere. 
As clearly shown in Fig. [TUJ the VACF of a microsphere 
in water is very different from exponential decay because 
of the hydrodynamic memory effects of water. 

The thin solid lines are expected to be correct for t ~ 
t c , and the thick sold lines are expected to be correct 
for t ^> t c . The intermediate regime t c < t < 100t c is 
poorly understood. A recent experiment has measured 
the VACF of a Brownian particle in water at VACF< 
0.35 0, [30|. A measurement of the VACF between 1 
and 0.35 is required in order to better understand the 
hydrodynamic effects and compressibility effects of water 
on Brownian motion 181-841. 



B. Experimental observation of the color of 
thermal force in a liquid 

Conventionally, the thermal force exerted on a Brown- 
ian particle are assumed to be a white noise. Due to the 
hydrodynamic memory effects of the liquid, this thermal 
force is in fact colored. Direct experimental observation 
of the color of the thermal force in a liquid was elusive 
until the recent work by Franosch et al [25|- They used 
optical tweezers with very large stiffness (k « 205 /iN/m 
for 2.9 fim diameter particles) to study the Brownian 
motion of a trapped microsphere. Their key insight is 
that at long times, strong trapping eventually dominates 
over friction and becomes the main force counteracting 
thermal force [27]. So the Langevin equation reduces to 
kx(t) ~ Ftherm(t)- Consequently, the correlations in the 
thermal force can be obtained by the position autocorre- 
lation function (F therrn (t)F therm (0)) « k 2 (x(t)x(0)). 

Figure fTTI shows the results of Franosch et aVs exper- 
iment on the color of thermal force [25|- The data in 
Fig. UTTa) clearly shows the departure of thermal force 
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from a white noise whose power spectrum is a horizon- 
tal line. The measured spectrum of the thermal force 
increases at higher frequencies. Franosch et al also ob- 
served resonances in Brownian motion in liquid where 
overdamped motion is often assumed. In order to ob- 
serve the resonance, they used optical tweezers with very 
large stiffness, and acetone as a liquid rather than water. 
The viscosity of acetone is about 3 times smaller than 
that of water. So the motion of a trapped bead is less 
damped in acetone than in water. The resonance can be 
clearly seen in Fig. Qjjb). Remarkably, this resonance is 
mainly due to the inertia of the liquid, rather than the 
inertia of the particle itself. 



(a) 1.7 




10 3 10 4 
Frequency, f (Hz) 

(b) L2 |L -r^^_ 1 




Frequency, f (Hz) 

FIG. 11. The color of thermal force (Figure adapted from Ref. 
[25l]). (a), The normalized power spectral density (PSD) of 
thermal force of an optically trapped melamine resin sphere 
(R = 1.45 /zm) in water (green circles) or acetone (blue cir- 
cles). The black lines are predictions of the hydrodynamic 
theory, (b), The normalized power spectral density of posi- 
tion. 



C. Experimental observation of the transition from 
ballistic to diffusive Brownian motion in a liquid 

Recently, Huang et al studied the Brownian motion 
of a single particle in an optical trap in water with sub- 
Angstrom resolution and measured the velocity autocor- 
relation function of the Brownian motion [24]. In their 
experiments, 77- due to the confinement of the optical 
tweezers was typically two orders of magnitude larger 
than Tf. So the role of the optical confinement can be 
neglected during the transition from ballistic to diffusive 
Brownian motion. 




10" 8 10 ' 1CT 6 io- 5 10 1 10 3 
Time (s) 



FIG. 12. Measured mean square displacement of silica spheres 
with diameters of 1 /am and 2.5 /xm (Figure adapted from Ref. 
[241]). The red lines with slope 2 show the expected behavior 
of ballistic Brownian motion. 




Time (s) 

FIG. 13. Normalized experimental VACF for a 2 /im di- 
ameter resin sphere in an optical trap (Figure adapted from 
Ref. [U]). The blue line shows the fitting of the Clercx-and- 
Schram theory (Eq. ([26]) ) to the whole experimental VACF. 
The grey (dashed) and green lines show the expectations of 
VACF for the same particle if the inertia of the liquid is ig- 
nored, in the absence (grey) and presence (green) of the har- 
monic trap. 
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Figure [12] shows Huang et aVs results of MSD's for 
lam and 2.5 /im silica particles from 10 ns to 1 ms 
|24j . The MSD's increase at short time scales, and reach 
plateaus at long time scales due to the confinement of the 
optical trap. At short time scales, the particle undergoes 
free, but correlated Brownian motion because of the iner- 
tia of the particle and the surrounding liquid. The effect 
of the particles' inertia on Brownian motion is visible in 
Fig. [121 The slopes of the MSDs are close to 2 at short 
time scales, as predicted for ballistic Brownian motion. 
It is remarkable that Huang et al achieved a temporal 
resolution of about 10 ns in the case of a 1 /im silica par- 
ticle. At this temporal resolution, Huang et al resolved 
a MSD as small as 0.0005 nm 2 , corresponding to 20 pm 
spatial resolution. The temporal resolution shown in Fig. 
[12] is about 0.1 r p for a 1 jam particle, and about 0.15 r p 
for a 2.5 /im particle. 

Figure [13] shows Huang et aUs experimental results of 
the VACF of a 2 /im resin particle and a theoretical 
fit with the Clercx-and-Schram theory (Eq. ((26])) |zi. 
The predictions of VACF from Langevin equations ne- 
glecting the hydrodynamic effects of water are also dis- 
played for comparison. The experimental results agree 
with the Clercx-and-Schram theory. For times shorter 
than Tf = 1 /is, the velocity correlations are smaller than 
the predictions of the Langevin equation neglecting the 
inertia of the fluid. At longer times, the correlations are 
stronger because of the vortex developed in the fluid. The 
small anti-correlation dip near 300 /is is because of the 
optical trap. 

Although the spatial resolution of Huang et aVs exper- 
iment is sufficient to observe the transition from ballistic 
to diffusive Brownian motion in MSD and even compute 
a VACF, it is not enough to measure the instantaneous 
velocity of Brownian motion in a liquid. The reason is 
that both MSD and VACF are insensitive to white noise 
of the detection system since they are averages over large 
data sets. On the other hand, the measurement of the 
instantaneous velocity is very sensitive to the detection 
noise. We will discuss about this issue in the following 
section. 



IV. EFFECTS OF DETECTION NOISE ON 
STUDYING BROWNIAN MOTION 

In the presence of detection noise [2(J HH, [24| , the mea- 
sured position of the microsphere can be expressed as 

(t) = x p (t)+x n (t), (31) 

where x p (t) is the real position of the microsphere, and 
x n (t) is the noise of the detection system. The mean 
square displacement (MSD) of the measured positions is 




MSD m8r (t) ([x msr (t -ft) 

(*o)] 2 > (32) 
= ([x p (t -ft)- x p (t )} 2 ) + ([x n (t -ft) - x n (t )} 2 ) 
+2{[x p (t -ft) - x p (t )] • [x n (t + 1) - x n (t )}) 
= MSD p (t) + MSD n (t), 

where the average is taken over all possible to- This 
derivation assumes no correlation between the real po- 
sition of the microsphere and the detection noise. In this 
case, the real MSD of the microsphere can be obtained 
by MSD p (t) = MSD^M) - MSD n (t), as is done by 
Huang et al in Ref. [18|, [24[. MSD n (t) is usually inde- 
pendent of time, as shown in Fig. [4] 

The measured velocity of the microsphere is 

r x x msr(t + ~y) ~ %msr(t ^r) , oq x 
Vmsr(t) ^ ) 

X p (t -f ~2~) X p (t ~2~) «^n(t "f ~2~) (t ~2~) 

= At + At 

= v p (t) +v n (i), 

where At <C r p . Because the measured velocity con- 
tains a noise signal v n (i), the smallest At of the detec- 
tion system may not be the best value for measuring 
the velocity. The data acquisition (DAQ) card creates 
noise when it converts an analog signal to a digital signal 
due to the finite number of bits. The minimum value 
of x n (t -f ^) — x n (t — ^) is limited by the DAQ card, 
thus v n (t) may be larger than the real velocity of the 
microsphere (v p (t)) if At is too small. 

The measured velocity represents the real instanta- 
neous velocity of the microsphere if At <C r p and v n (t) is 
negligible. This requires (t^ sr ) ^> (v^)- One can check 
whether this condition is satisfied by comparing the sig- 
nal when a microsphere is trapped in the optical tweezer 
and when there is no microsphere in the optical tweezer. 
The relation between (^ sr ) and (t> 2 ) can be obtained 
from the measured MSD's: 

i 2 \ _ i [ x msr(t + — ) — X msr (t 2^)] 2 > , 
\ v msr) \ At 2 

_ MSD p (At) MSD n (At) 
~ At 2 + At 2 

Thus (v^) > (vl) is equivalent to MSD msr (At) > 
MSD n (At). 

The measured velocity autocorrelation function is 

(Vmsr(t + t )v m sr(to)) = (v p (t + t )v p (t )) (35) 

-f (v n (t + t )v n (to)). 

If the noise of the detection system has almost no cor- 
relation (white noise), the last term of this equation can 
be neglected. Thus 

(Vmsr(t + to)v m8 r(to)) = {v p (t + t )v p (t )}. (36) 

So the measurement of the velocity autocorrelation func- 
tion is not sensitive to the noise of the detection system 
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[24J . On the other hand, the measurement of the in- 
stantaneous velocity is very sensitive to the noise of the 
detection system [21] . 

If the detection system samples the position of the mi- 
crosphere every dt that is much shorter than the required 
temporal resolution At, we can reduce the noise in the 
measured velocity by using successively averaged posi- 
tions to calculate the velocity. Let At = N dt, where 
N < Tp/dt, then 



1 N 



~(t + jdt 



7V + 1 



dt). (37) 



The measured velocity becomes 



-(*) 



Then the velocity noise is 



N 



At) 
2 J 



At 



(38) 



W = WnH E + ^ ~ + ^) ( 39 ) 



iV 



-2_^Xn{t + jdt — dt - — )]. 

On average, the rms amplitude of f n (£) in Eq. ([39]) is 
N\/N times smaller than that of [x n (£ + dt/2) — x n (t — 
dt/2)]/dt if the position noise x n (t) is a white noise. If 
N = 10, the noise can be reduced by a factor of 32, which 
is a significant number. The condition of using the av- 
eraging method is that dt <C r p /N. So the detection 
system must have a very high temporal resolution. This 
averaging method has been utilized in both the measure- 
ment of the VACF of a Brownian particle in a liquid [24| 
and the measurement of the instantaneous velocity of a 
Brownian particle in air [21]. 



V. FUTURE 

As shown in this review, optical tweezers have be- 
come an indispensable tool for studying Brownian mo- 
tion at short time scales. The instantaneous velocity of 
a Brownian particle trapped in a gas has been measured 
[21| . Recently, the velocity autocorrelation function of a 
Brownian particle in water was measured successfully for 
M*M2)>< 0.35 |H MM- 

The instantaneous velocity of a Brownian particle in 
water is much more difficult to measure, and has not 
been measured to date. A successful measurement of 
the instantaneous velocity of a Brownian particle in a 
liquid will complete the task that was considered by Ein- 
stein more than 100 years ago [ll|, [l2| and open a new 
door for studying Brownian motion. For example, the 
results can be used to test the modified energy equipar- 
tition theorem (v(0)v(0)} = k B T '/M* and the Maxwell- 
Boltzmann velocity distribution. A measurement of the 
VACF > 0.35 will deepen our understanding of the hy- 
drodynamic effects and compressibility effects of a liquid 
on Brownian motion [8ll-[84) . The ability to measure the 
instantaneous velocity of a Brownian particle will be in- 
valuable in studying nonequilibrium statistical mechanics 
[85l |86| . The Brownian motion of a suspended particle 
can be used for microrheology to probe the properties of 
fluids, such as viscoelastic fluids [87-89], and surrounding 
environments |9CH92j. 
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